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Abstract 

C 4 plants have a biochemical carbon-concentrating mechanism that increases C0 2 concentration around Rubisco in the 
bundle sheath. Under low light, the activity of the carbon-concentrating mechanism generally decreases, associated 
with an increase in leakiness ((])), the ratio of C0 2 retrodiffusing from the bundle sheath relative to C 4 carboxylation. This 
increase in <\> had been theoretically associated with a decrease in biochemical operating efficiency (expressed as ATP 
cost of gross assimilation, ATP/GA) under low light and, because a proportion of canopy photosynthesis is carried out 
by shaded leaves, potential productivity losses at field scale. Maize plants were grown under light regimes representing 
the cycle that leaves undergo in the canopy, whereby younger leaves initially developed under high light and were then 
re-acclimated to low light (600 to 1 00 uE-m" 2 s" 1 photosynthetically active radiation) for 3 weeks. Following re-acclimation, 
leaves reduced rates of light-respiration and reached a status of lower q>, effectively optimizing the limited ATP resources 
available under low photosynthetically active radiation. Direct estimates of respiration in the light, and ATP production 
rate, allowed an empirical estimate of ATP production rate relative to gross assimilation to be derived. These values were 
compared to modelled ATP/GA which was predicted using leakiness as the sole proxy for ATP/GA, and, using a novel 
comprehensive biochemical model, showing that irrespective of whether leaves are acclimated to very low or high light 
intensity, the biochemical efficiency of the C 4 cycle does not decrease at low photosynthetically active radiation. 

Keywords: Bundle sheath, A 13 C, irradiance, isotopic discrimination, leakiness, low light, mesophyll, efficiency, PPFD. 



Introduction 

The C 4 pathway of photosynthesis has been attracting increas- 
ing interest in recent years for its high crop productivity potential 
in the face of global warming and population pressure (Friso 
et aL, 2010; Zhu et aL, 2010; Covshoff and Hibberd, 2012). C 4 
photosynthesis evolved from C 3 photosynthesis under the envi- 
ronmental pressure of declining ambient C0 2 and increasing 
transpiration demand in semi-arid environments (Griffiths et aL , 
2013; Osborne and Sack, 2012). Under optimal conditions, 



characterized by high temperatures and high light intensi- 
ties, C 4 plants have higher photosynthetic rates than C 3 plants 
(Ehleringer and Pearcy, 1983; Pearcy and Ehleringer, 1984) and 
very high productivity. Many C 4 plants have been domesticated 
and represent irreplaceable sources of food, biomass, and bio- 
energy. For instance, maize (Zea mays, L.), a C 4 plant of the 
NADP-malic enzyme (NADP-ME) subtype, is the leading grain 
production cereal (www.fao.org/statistics). 



Abbreviations: BS, bundle sheath; CCM, carbon-concentrating mechanism; ETR, electron transport rate; HL, plants grown under high light; HLLL, mature leaves 
re-acclimated under low light; IRGA, infrared gas analyser; LCR light compensation point; LL, plants grown under low light; PAR, photosynthetically active radiation; 
PEP, phosphoenolpyruvate; PEPC, phosphoenolpyruvate carboxylase; PEPCK, phosphoenolpyruvate carboxykinase; PGA, 3-phosphoglyceric acid; PPDK, pyru- 
vate phosphate dikinase; PSI, photosystem I; PSII, photosystem II; RGB, red, green, and blue; RuBR ribulose 1 ,5-bisphosphate. 

© The Author 2014. Published by Oxford University Press on behalf of the Society for Experimental Biology. 

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.Org/licenses/by/3.0/), which 
permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. 
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The high productivity of C 4 plants results from ana- 
tomical and biochemical differentiation of the leaf paren- 
chyma. Externally mesophyll cells and internally bundle 
sheath (BS) cells are coupled to operate a biochemical 
carbon-concentrating mechanism (CCM) that increases 
the C0 2 concentration in BS, the cellular compartment 
where Rubisco is exclusively expressed, resulting in active 
suppression of the oxygenase activity of Rubisco. Since 
BS and mesophyll cells are connected by plasmodesmata, 
some C0 2 retrodiffuses (C0 2 leakage). The extent of 
C0 2 retrodiffusion is still debated, but it is accepted that 
the permeability to C0 2 diffusion (BS conductance, g BS ; 
Table 1) varies between different species and individual 
plants. A useful term to describe this concept, which was 
coined by Farquhar in the description of carbon isotope 
discrimination (Farquhar, 1983) is leakiness (</>), defined 
as the rate of C0 2 retrodiffusing (leak rate) relative to the 
phosphoenolpyruvate (PEP) carboxylation rate (V P ). Since 
Rubisco C0 2 fixation (in BS) is complementary to leakage 
(out of BS), 4> can be used as a proxy for the coordination 
between the CCM and C 3 assimilatory activity (Henderson 
et al., 1992; von Caemmerer, 2000; Tazoe et al, 2006, 2008; 
Kromdijk et al, 2010; Ubierna et al, 2011; Bellasio and 
Griffiths, 2013). 

The CCM has a notable metabolic cost: out of the the- 
oretical minimum of five ATP molecules required for the 
gross assimilation of one C0 2 , two ATPs are consumed by 
the CCM (Furbank et al, 1990; Bellasio and Griffiths, 2014) 
in the costly regeneration of PEP. The common interpreta- 
tion of C 4 physiology assumes that, at steady state, leaking 
C0 2 is entirely refixed by PEP carboxylase (PEPC); hence, 
anatomical features are tightly bound to biochemical and 
energy traits. Plants with a higher g BS would have higher rate 
of C0 2 retrodiffusion, increased CCM cost, and a higher 
ATP demand for gross assimilation (ATP/ GA), which is the 
overall biochemical operating efficiency of C 4 photosynthe- 
sis. For this reason, 0 has been used to derive ATP/GA (von 
Caemmerer, 2000; Tazoe et al, 2008); for instance, plants 
with higher 4> are considered to have higher ATP/G^4 and, 
therefore, lower biochemical operating efficiency. We will 
show that these assumptions hold true only under high light 
intensities. 

Because of these anatomical, biochemical, and energetic 
complexities, C 4 metabolism is highly sensitive to limit- 
ing light intensities (see Ubierna et al., 2011 for review). 
Recently, studies have focused on characterizing the pro- 
gressive increase in carbon isotope discrimination that is 
usually seen as light intensity decreases at both leaf (Tazoe 
et al, 2008; Kromdijk et al, 2010; Pengelly et al., 2010; 
Bellasio and Griffiths, 2013; Ubierna et al, 2013) and can- 
opy (Kromdijk et al, 2008) levels. The theoretical consid- 
erations highlighted above have associated this increase in 
4> with decreased C 4 efficiency and a potential loss of pho- 
tosynthetic carbon uptake (Furbank et al., 1990; Kromdijk 
et al., 2008; Tazoe et al, 2008). Empirical evidence was 
needed to validate this suggestion and to explore the strate- 
gies that mature C 4 leaves deploy to cope with reduced light 
intensities. 



Low light responses are highly relevant for C 4 canopy pro- 
ductivity, since up to 50% of net C0 2 uptake (Baker et al, 
1988; Long, 1993) is fixed by shaded leaves, under a light 
intensity which is typically one-twentieth of full sunlight 
(Shirley, 1929). In a forest canopy leaves are subjected to a 
similar degree of exposure throughout the year, whereas 
in crop canopies most fully expanded leaves progressively 
acclimate to shade under newly emerging leaves. This long- 
term acclimation is accompanied by transitory, short-term 
responses such as daily shading, or more transient sunflecks. 
Furthermore, there is a gradient of leaf age down the crop 
canopy, with younger leaves exposed to full sunlight at the 
top of the canopy and older leaves subsequently exposed to 
canopy-filtered light. 

Previously we studied how long-term acclimation to 
low light influenced short-term responses to illumination 
(Bellasio and Griffiths, 2013). Plants grown under low light 
(LL) showed a capacity for maintaining low 4> even under 
decreasing light intensities, whereas <p increased in equiv- 
alent plants grown under high light (HL). We suggested 
several mechanisms whereby C 4 leaves adapted through- 
out growth to low-light conditions could maintain high 
photosynthetic conversion efficiency during steady-state 
photosynthesis. 

In this study we grew maize plants under a light regime 
representing the acclimation of leaves shaded by an over- 
growing canopy, consisting of 3 weeks under high light 
followed by 3 weeks under diffuse, low light intensity. The 
leaf-level ATP-production rate (J A tp) was derived from gas 
exchange measurements under low 0 2 in combination with 
photosystem II (PSII) photochemical yield, measured C0 2 
assimilation rate, and online isotopic discrimination during 
photosynthesis (A). A full isotopic discrimination model was 
used to derive cf> from A (Farquhar, 1983; Ubierna et al. , 20 1 1 ; 
Farquhar and Cernusak, 2012). With the directly derived val- 
ues for J ATP , the empirical ATP cost of gross and net assimi- 
lation (respectively, J ATP IGA and J ATP IA) could be calculated 
and compared with the predicted ATP cost of assimilation 
(ATP/G^4). Mature leaves that had re-acclimated under low 
light (HLLL) showed very similar traits to LL plants. HLLL 
plants deployed two strategies to optimize the scarce ATP 
resources under low light: (i) the reduction of respiration 
in the light (R U ght) and (ii) the reduction of leakiness (0). 
The comparison of J ATP IGA with ATP/GA estimated with a 
novel metabolic model showed that C 4 photosynthetic effi- 
ciency was constant in the vicinity of the light compensation 
point (LCP): thus, the predicted decrease in biochemical con- 
version efficiency based on (p increasing under limiting light 
does not occur. 



Materials and methods 

Plants 

Plants were grown at the Plant Growth Facility located at the 
University of Cambridge Botanic Garden in controlled-environ- 
ment growth rooms (Conviron, Winnipeg, Canada) set at 16-h 
day length, 25/23 °C (day/night), and 40% relative humidity. The 
growth protocol was designed to standardize age and watering 



Table 1 . Definitions, equations, and variables used 
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Symbol 



Definition 



Values/units/references 



ATP/G4 

b 3 



b 3 
b< 



c, 

Cm 



E 

F 

f 

GA 



9bs 
9t,< 
9s 
Jatp 

J ATP /A 

J ATP /GA 

Jmod 



Net assimilation 

13 C fractionation due to diffusion of C0 2 in air. Due to vigorous 
ventilation we ignored fractionation at the boundary layer. 
13 C fractionation due to diffusion of C0 2 in water 
Predicted ATP demand for gross assimilation, i.e. predicted 
biochemical operating efficiency 

13 C fractionation during carboxylation by Rubisco including 

e'-Rimw + f-F 



respiration and photorespiration fractionation b 3 =b 3 '- 



V. 



13 C fractionation during carboxylation by Rubisco 
(excluding respiration and photorespiration fractionation) 
Net fractionation by C0 2 dissolution, hydration and PEPC 

, s'Rm 
V P 



carboxylation including respiratory fractionation b 4 = b 4 



Net fractionation by C0 2 dissolution, hydration, and PEPC carboxylation 
(excluding respiratory fractionation) 



xJ 



ATP n UGHT 



CO, concentration in the BS; C F 



- + Cm 



9bs 



C0 2 concentration in the intercellular spaces as calculated by the IRGA. 

A 



C0 2 concentration in the mesophyll; C M = C r 



13 C fractionation during decarboxylation 



9m 



3 C fractionation during decarboxylation, including the correction for 



measurement artefacts: e =e + <5 13 C, 



measurements 



-5 ^ 3 C gr0Wthcnl}mt)er although 



there may be some error at low light intensities if recent photosynthate is 

not the substrate for respiration 

13 C fractionation during internal C0 2 dissolution 

transpiration rate (calculated by the IRGA software, parameter Trmmol) 
Rate of photorespiratory C0 2 evolution F = 0.5 ■ V 0 

13 C fractionation during photorespiration 

Gross assimilation GA = A + R UGHT 

conductance to diffusion of C0 2 in air (calculated by the 

IRGA software, parameter CndC02) 

BS conductance to C0 2 , calculated by fitting J M0D to J ATP 

Mesophyll conductance to C0 2 

Stomatal conductance to C0 2 



ATP production rate J A 



3 GA 



Y(ll) 
0.59 Y(ll) Low02 



Low Oi 



ATP production rate relative to net assimilation 
ATP production rate relative to gross assimilation 



Modelled ATP production rate J mod 

9bs 



-y + Vy 



-Awz 



2w 



y 



1-x 
3 



A 



+ Cm 



Rm 

9bs 



ay 



0.047 



1 + 



where: w 
Rught 



x-x 

6A 



1 + - 



_ r 7 9bsY C m , 

" 9BS^M 7, + (r>LIGHT 



(Rught+AMI 



lay 
' 3-0.047 



|imol-rrr 2 -s~ 1 

4.4%o (Craig, 1953; Kromdijk etal., 2010) 
0.7%o (O'Leary, 1984) 



% 0 (Farquhar, 1983; Ubierna et al. , 2013) 

30%o (Roeske and Oleary, 1984) 

% 0 (Farquhar, 1983; Henderson et al. , 1992) 



-5.7% 0 at 25 °C but variable with temperature 
(Farquhar, 1983; Henderson etal., 
1992; Kromdijk etal., 2010) 
nmol-mol -1 



Hmol-mol -1 (Li-cor 6400 manual eqn 1.18) 
nmol-mor 1 

Oto -10%o (Gillon and Griffiths, 1997; Ghashghaie 
ef al. , 2001 ; Igamberdiev ef al. , 2004; Hymus et al. , 
2005; Barbour etal., 2007; Sun et al., 2012); -6% in 
this study (Kromdijk ef al. , 201 0) 



% 6 13 C, 



measurements " 



-6.38%o; 6 13 C QI 



3%o 



(Wingate ef al. , 2007) 



I . 1 %» (Vogel ef al. , 1 970; Mook ef al. , 1 974; Vogel, 
1980) 

mmol-rrf 2 -s~ 1 

|xmol-m~ 2 -s _1 (von Caemmerer, 2013; N. Ubierna, 
personal communication) 

I I . 6%o (Lanigan ef al. , 2008) 
|imol-rrr 2 -s~ 1 
mol-m _2 -s _1 

mol-rrr 2 -s- 1 (Bellasio and Griffiths, 2013) 
1 mol-rrr 2 -s- 1 -bar 1 (Kromdijk ef al., 2010) 
mol-m" 2 -s" 1 

|imol-m _2 -s" 1 (Bellasio and Griffiths, 2013) 



ATP / C0 2 
ATP / C0 2 

|imol-m~ 2 -s _1 (von Caemmerer, 2000; Bellasio and 
Griffiths, 2013; Ubierna etal., 2013) 
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Symbol 



Definition 



Values/un its/references 



0 BS 
O,, 

RuGHT 

Rm 

s 
f 



0 2 mol fraction in the BS cells (in air at equilibrium) O ss = 0 M + 



aA 
0.047 g BS 



0 2 mol fraction in the mesophyll cells (in air at equilibrium) 
Respiration in the light 

Mesophyll non-photorespiratory C0 2 production in the light R M = 0.5 R U qht 

Fractionation during leakage of C0 2 out of the BS cells 
(1 + a) E 

Ternary effects f = 



2000 g m 



Rubisco carboxylation rate V c 



Rubisco oxygenation rate V 0 



= (A + R U cht) 
1 y'Obs 

0.5 



^mol-mor 1 (von Caemmerer, 2000) 

210000 nmol-mol" 1 
nmol-rrf 2 -s~ 1 

nmol-m" 2 -s" 1 (von Caemmerer, 2000; 
Kromdijk ef a/. , 201 0; Ubierna ef a/. , 201 3) 
1 .8%o (Henderson ef a/. , 1 992) 
%o (Farquhar and Cernusak, 2012) 

nmol-rrr 2 -s~ 1 (Ubierna era/., 2011) 



nmol-m 2 -s 1 (Ubierna ef a/. , 2011) 



PEP carboxylation rate V P 



Y(ll) 



J ATP partitioning factor between C 4 activity (V P ) and C 3 activity V c +V a 
(reductive pentose phosphate pathway and photorespiratory cycle) 



Yield of photosystem II Y(ll) 



Set at 0.4 (von Caemmerer, 2000; Kromdijk ef a/., 
201 0; Ubierna ef a/. , 201 1 , 201 3), except for the 
calculation of eqn 2 where x was not constrained, 
dimensionless (Genty ef a/., 1989) 



V 
A 



6 13 C 



Fraction of PSII active in BS cells 

Half of the reciprocal of the Rubisco specificity 
13 C Isotopic discrimination A - 0 ^ 



1 + <5„ 



where: % - 



C fl — C, 



see supporting 



Fig. S1 ; 6 e is the isotopic composition of the reference gas. 5 0 is the isotopic 
composition of the gas leaving the cuvette. C e and C 0 represent the C0 2 mole 
fraction respectively entering and leaving the cuvette corrected for differing 
amounts of water vapour according to (von Caemmerer and Farquhar, 1 981 ). 

13 C isotopic composition relative to Pee Dee Belemnite 
Leakiness; defined as the leak rate relative to V P It was estimated 
with the isotope method including respiratory and photorespira- 
tory fractionation, ternary effects and estimating C BS with the C 4 model 
C BS -C M b 4 C M (1 + f) + a(C a -C i )-C a A 0SS (1-f) 

Cm (1 + f)[C a A 0SS (1-f)-a(C a -C^-b.Cgs +s(C BS -C M )\ 



0.15 (von Caemmerer, 2000; Edwards and Baker, 
1993; Kromdijk ef a/. , 2010) 
0.000193 (von Caemmerer, 2000) 
%o (Evans ef a/., 1986) 



Too 

dimensionless (Farquhar and Cernusak, 201 2) 



conditions throughout the experiment. Two light environments 
were established — high-intensity direct light (photosynthetically 
active radiation, PAR = 600 u,E-irr 2 -s _1 ) and low-intensity diffuse 
light (PAR = 100 u.E-m -2 -s ) — obtained using shading to mimic 
the understory of a canopy. Maize seeds (Zea mays L. ¥ l hybrid 
PR31N27; Pioneer Hi-bred, Cremona, Italy) were sown weekly in 
1.5 1 pots filled with Levington pro M3 pot and bedding compost 
(Scotts, Godalming, Surrey, UK). Plants were grown in three sets of 
conditions: (i) HL plants were grown for 3 weeks under high light 
(fully expanded fourth leaf stage); (ii) LL plants were grown for 4 
weeks under low light (fully expanded fourth leaf stage); and (iii) 
HLLL plants were grown for 3 weeks under high light, the youngest 
fully expanded leaf was marked, and then plants were grown for the 
following 3 weeks under low light. Plants were manually watered 
daily, with particular care to avoid overwatering. When ready, 
plants were measured once and then discarded. Measurements 



were performed on the youngest fully expanded leaf of HL and LL 
plants, and on marked leaves of HLLL plants. 

Gas exchange measurements with concurrent PSI/PSII yield 
and online carbon isotopic discrimination (A) 

The experimental setup was previously described in detail (Bellasio 
and Griffiths, 2013). Briefly, an infrared gas analyser (IRGA; an 
LI6400XT, Li-cor, Lincoln, NE, USA) was fitted with a 6400-06 
PAM2000 adapter and with a Li-cor 6400-18 red, green, and blue 
(RGB) light source. RGB light was used because, by providing equal 
fractions of R, G, and B, it is likely to distribute excitation between 
mesophyll and BS cells with a more similar pattern to natural white 
light than the conventional 90% R/10% B source. The IRGA was fed 
with C0 2 (6 13 C = -6%o; Isi Soda, Vienna, Austria) and either a mix- 
ture of 2% Q 2 /N 2 or ambient air. Photosystem I (PSI) yield and PSII 
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yield ( Y( II); see Table 1) were measured using a Dual Pam-F (Heinz 
Walz GmbH, Effeltrich, Germany). Pulse intensity was set to 20 
mE*m~ 2, s , enough to saturate fluorescence and PSI signals (which 
occurred between 8 and 10 mE-m~ 2 -s -1 ; data not shown). The block 
temperature was set at 26 °C so as to maintain the leaf temperature 
close to 25 °C. The IRGA was connected to a cryogenic H 2 0- and 
C0 2 -trapping purification line. Each day, one plant was subject to a 
RGB-light-response curve, under 2% 0 2 and C a = 600 (xmol-mol -1 
[to determine the relationship between electron transport rate (ETR) 
and J A tp] an d a second RGB-light-response curve under 21% 0 2 and 
reference C0 2 set at 400 urnd-moF 1 , during which exhaust gas was 
trapped to determine A. With this procedure each day the 6 13 C com- 
position of a total of 12 C0 2 samples and six C0 2 references (repre- 
senting responses to decreasing irradiances of one individual plant) 
were analysed directly using a VG SIRA dual-inlet isotope ratio 
mass spectrometer (modified and maintained by Pro-Vac Services, 
Crewe, UK). A was calculated as reported in Table 1 (Evans et al, 
1986). Y( II) was determined at each light level for both light curves. 
J ATP was calculated individually at each irradiance by multiplying 
the relationship between ETR and J ATP (determined at low 0 2 ) by 
the ratio between Y(II) at ambient and low 0 2 (Table 1). R light \n&& 

calculated as the j'-intercept of the linear regression of net assimila- 
tion, A, against PAR . (Table 1; Yin et al, 2011b; Bellasio 

and Griffiths, 2013). Although we did not find significant differ- 
ences with values for dark respiration (measured with the IRGA 
every 10 s for 4min and averaged, with flow rate of 50 umol-s" 1 ) 
or with values of Rlight estimated through non-linear curve fit- 
ting of light-response curves (non-rectangular hyperbola; Prioul 
and Chartier, 1977; Dougherty et ah, 1994), we preferred the linear 
curve fitting described above for the robustness and the simplicity 
detailed in Yin et al (2011a). The LCP was calculated using dedi- 
cated software (Photosyn assistant 1 .2, Dundee Scientific, Dundee, 
UK). Assimilatory light-response curves were transformed logarith- 
mically and subject to analysis of variance (ANOVA; Genstat). The 
transformation was necessary to normalize the residuals and thereby 
avoid the artefactual intepretation of significance (i.e. significant 
differences only at higher light intensities). Responses to decreas- 
ing light intensities were subject to repeated-measures ANOVA 
(Genstat); point estimates were subject to ANOVA and Tukey mul- 
tiple comparisons as appropriate (Genstat). 

Leakiness cf> from isotopic discrimination A 

Modelling was previously described in detail (Bellasio and Griffiths, 
2013), and equations are reported in Table 1. Briefly, leakiness, <p 
was resolved from A using the full model of Farquhar, as recently 
integrated to take into account the 'ternary' effects, i.e. the effect of 
water molecules diffusing outward stomata on C0 2 molecules dif- 
fusing inwards through still air (Farquhar and Cernusak, 2012). In 
this model, the weighted individual fractionations of the discrimi- 
nating processes operating in C 4 photosynthesis are summed. This 
model requires the C0 2 concentration in the different cellular com- 
partments (notably mesophyll and BS cells), which were calculated 
by means of the validated C 4 photosynthesis model, in the light- 
limited form (later 'C 4 model'; von Caemmerer, 2000). The C 4 model 
was in turn parameterized with the light-response data (A, Q, C m 
J AT p) and Rlight- BS conductance, required to parameterize the C 4 
model, cannot be measured directly but it can be estimated by fitting 
the C 4 model to a measured quantity. In the A/A' fitting (Kromdijk 
et al, 2010; Ubierna et al, 2013), the C 4 model is rearranged to 
express a modelled isotopic discrimination and fitted to values for 
A. Here, we used the 'J/J' fitting, which we have recently described 
(Bellasio and Griffiths, 2013), whereby the C 4 model is rearranged 
to express a modelled ATP production rate J mod ar >d fitted to the 
empirically derived estimate for the leaf-level ATP-production rate 
J AT p, described above. This procedure yielded a value for g BS for 
each individual plant which was obtained independently from A, 



and did not suffer the circularity of the A/A' fitting, arising from 
calculating g BS and leakiness from the same values for A (Bellasio 
and Griffiths, 2013). 

Empirical and predicted ATP cost of gross assimilation 

We refer to empirical ATP cost of net and gross assimilation as 
J A tpIA and J ATP /GA, while we refer to predicted ATP cost of gross 
assimilation as ATP/ GA. 

The empirical ATP cost of net and gross assimilation was calcu- 
lated from the data obtained during the experiment. Firstly, the meas- 
ured leaf-level ATP cost of net assimilation (J ATP IA) was calculated 
from J ATP and net assimilation, A. The derivation of J ATP through the 
low 0 2 -ETR method was described above (see also Table 1). J AP plA is 
relevant to net productivity and shows how much ATP the plant has 
to spend for net gain of a C0 2 molecule. Then, the leaf-level ATP cost 
of gross assimilation (J ATP IGA) was calculated using values for GA, 
derived by summing A plus Rlight calculated by curve fitting (see 
above). J ATP IGA is relevant to C 4 biochemistry and shows the empiri- 
cal conversion efficiency of C0 2 into sugars. It is worth stressing that 
these values for J ATP IGA are derived with a novel method based on 
gas exchange under low 0 2 (Yin et al, 2011a, 2011b; Bellasio and 
Griffiths, 2013). The low 0 2 -ETR method relies on two assumptions. 
First is that the partitioning of NADPH to photosynthesis does not 
change between ambient and low 0 2 . This is a fair assumption since 
NADPH use by alternative sinks (e.g. nitrogen reduction) is generally 
dependent on light intensity and hence it is only marginally influ- 
enced by 0 2 partial pressure (Yin et al. , 2004, 2009; Yin and Struik, 
2012). Second, Rlight does not vary between low and ambient 0 2 . 
This is also a fair assumption because any 0 2 effect is generally neg- 
ligible (Badger, 1985; Gupta et al, 2009). The low 0 2 -ETR method 
does not rely on the assumptions used in the traditional derivation 
based on leaf absorptance and PSII optical section (von Caemmerer, 
2000) and should therefore better represent the actual biochemical 
ATP demand of the portion of leaf subject to ecophysiological char- 
acterization. Because of the difficulty in deriving J ATP IGA based on 
leaf absorptance, and the difficulty in capturing the stoichiometry at 
the electron transport chain, the ATP cost of gross assimilation has 
often been predicted (e.g. Tazoe et al , 2008). 

A traditional way to predict ATP/ GA uses leakiness as the sole 
proxy (4> approach; Furbank et al, 1990; von Caemmerer, 2000; 
Tazoe et al, 2008). The (j) approach relies on the assumption that the 
ATP cost of the C 3 activity is invariably 3 ATP/C0 2 (photorespira- 
tion is neglected) while the ATP cost of the CCM depends solely on 
(j>. This implies that the CCM is driven solely by the activity of PEPC 
and that all the retrodiffusing C0 2 is refixed. Under these assump- 
tions the ATP cost of the CCM is calculated by multiplying the over- 
all ATP cost of PEPC (2 ATP/C0 2 ) by the ratio of C0 2 overcycling 
[1/(1-0)]. The total ATP cost of gross assimilation results from sum- 
ming the cost of the C 3 activity plus the cost of the CCM (see eqn 5 
in Tazoe et al 2008, or eqn 4.55 in von Caemmerer, 2000): 

ATP „ 2 

=3 + (1) 

GA 0 1-0 v 

Here the subscript (j> recalls that ATP/GA is derived from leakiness. 
Eqn 1 was solved for the three types of plants (HL, LL, and HLLL) 
and light intensities from 50 to 500 nE-m~ 2 -s using the values of <p 
derived from isotopic discrimination. 

We propose a different approach to estimate ATP/ GA, whereby the 
ATP demand of all biochemical processes underpinning assimilation 
(hence B approach) are summed. The B approach is comprehensive, 
and requires the quantification of all processes contributing to C 4 
photosynthesis. We used the validated C 4 model (von Caemmerer, 
2000), as recently integrated to describe the C 4 energetics (Bellasio and 
Griffiths, 2014). The biochemical processes considered are: 3-phos- 
phoglyceric acid (PGA) reduction, starch synthesis, PEP regenera- 
tion, ribulose 1,5-bisphosphate (RuBP) regeneration, and glycolate 
recycling, while the PGA consumed by mitochondrial respiration is 
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subtracted as likely to be consumed by basal metabolism (for deriva- 
tion see Bellasio and Griffiths, 2014). ATP/GA B was calculated as: 

A TP 1 A 1 

— - = W C + -V 0 +- + PEPCK + 2PPDK--R LIGHT (2) 
KjA b 2o 3 

Where the subscript B recalls that all the biochemical processes were 
summed, V c is the Rubisco carboxylation rate, V Q is the Rubisco oxy- 
genation rate, A is net assimilation, PEPCK is the PEP carboxykinase 
(PEPCK) rate, and PPDK is the pyruvate phosphate dikinase (PPDK) 
rate. PEPCK was assumed to regenerate 20% of the PEP required by 
PEPC, which is close to the expected value under natural white light 
(Bellasio and Griffiths, 2014); the remainder was regenerated through 
PPDK. PEPC rate (V P ), V c , and V 0 were calculated with the validated 
von Caemmerer C 4 model (Table 1), in the light-limited form (von 
Caemmerer, 2000; Bellasio and Griffiths, 2013). The model was con- 
strained at each light intensity with the values for A and J ATP shown in 
Fig. 1, with the values for CjC a and C BS shown in Fig. 2, with the values 
for R L jght and gBs reported in Table 2, and with the values for <p shown 
in Fig. 3. A parameter, known as x, is required to solve the C 4 model 
(Table 1), which partitions the ATP available between the CCM activity 
and the C 3 activity (PGA reduction, RuBP regeneration, and glycolate 
recycling). For the purposes of these calculations, rather than using a 
fixed value of x, such as 0.4 (von Caemmerer, 2000; Kromdijk et al., 
2010), we allowed x to vary to get the best fit for the parameters above. 

Results 

Physiological response to decreasing light intensities 

Figure 1 shows the responses of maize plants to decreasing 
irradiance when grown under three different light regimes. 
Assimilation (^4) significantly differentiated plant responses 
(Fig. 1A). LL plants had the highest A at PAR lower than 
500 uE-m~ 2 -s _1 . HL plants had the highest A at saturating 
PAR and the lowest A at PAR lower than 250 LiE-m~ 2 -s _1 . 
HLLL plants had the lowest A under saturating light 
(although these leaves were now 3 weeks older), but as light 
decreased between PAR 250 and 0 uE-m~ 2 -s _1 the response 
approached that of LL plants. Consistently, the LCP and 
Rught of HLLL plants were similar to those of LL plants, 
and clearly lower than those of HL plants (Table 2). 

The total ATP production rate (J ATP ) is shown in Fig. IB. 
J ATP was derived from gross assimilation under low 0 2 and 
then corrected for photorespiration at ambient 0 2 using the 
ratio of photochemical yield. At high PAR, J ATP tracks the 
pattern of A; however, at low PAR, J ATP of all plants was sim- 
ilar, suggesting that the higher A of LL and HLLL plants at 
limiting PAR (inset in Fig. 1 A) was achieved through a higher 
conversion efficiency and lower respiration rate (Table 2). 
Isotopic discrimination during photosynthesis (A) is shown in 
Fig. 1C. In HL and HLLL plants A increased substantially at 
PAR lower than 250 LiE-nr 2 -s _1 , although in HLLL plants A 
was, on average, lower than for HL plants. LL plants showed 
a more gradual increase under decreasing PAR. 

Figure 2A shows stomatal conductance (g s ) and Fig. 2B 
shows CJC CI . CJC U differentiated clearly between growth condi- 
tions, and was lowest in LL plants and highest in HL plants, 
while HLLL plants had intermediate values at all levels of 
PAR. C/C fl was higher than 0.5 at PAR <125 LiE-nTV (LL 
plants), reflecting the efforts made during the measurements to 
induce stomatal opening. A high C/C„ was functional in the 
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Fig. 1. Maize responses to decreasing light intensities for plants grown 
under high light (HL), low light (LL) or LL following HL (HLLL). (A) Net 
assimilation (4). The curves were fitted to calculate the LCP (Table 2). The 
inset shows a magnification at the lowest PAR. (B) Total ATP production rate 
(Jatp)> measured with the low 0 2 -ETR method (see Materials and methods 
section on gas exchange measurements). (C) Online isotopic discrimination 
during photosynthesis (A). Error bars represent one SE (n = 6). 

resolution of the isotopic discrimination model, to maximize 
the contribution of biochemical processes over the stomatal 
contribution to total isotopic discrimination (Table 1; Cernusak 
et a!., 2013). This was especially important for HL plants which 
have, under low light, lower assimilation than LL plants (and 
higher |, Table 1, Fig. SI; Evans et al., 1986). Figure 2C shows 
the C0 2 concentration in BS (C BS ), which was estimated by fit- 
ting a C 4 photosynthesis model, rearranged to express J mod, to 
the values for J ATP described above. The difference between con- 
ditions was not significant, and was due to a small difference in 
permeability to C0 2 retrodiffusion out of the BS (g BS ; Table 2). 

Leakiness 

Figure 3 shows leakiness, 4>, over the experimental range of 
PAR. These values were resolved from A through a full iso- 
topic discrimination model (Farquhar and Cernusak, 2012), 
parameterized using the C 4 model and fitted to J ATP , using 
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the recently described J/J fitting (Bellasio and Griffiths, 
2013). cp significantly differentiated the three types of plants 
(P = 0.008). HL plants had higher cp than LL and HLLL 
under limiting PAR, with <p increasing from 0.25 to 0.35 
under decreasing PAR. HLLL plants had the lowest <p under 
light intensities higher than 75 uF-m~ 2 -s _1 . Under low light 
intensities they showed an increase in <p with a similar trend 
to that of HL plants. In LL plants (p was close to 0.24 and 
only marginally affected by light intensity. 
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Fig. 2. (A) Stomatal conductance and (B) C/C a responses to decreasing 
light intensity, under different light qualities, for plants grown under high light 
(HL), low light (LL), or LL following HL (HLLL) measured by gas exchange. 
(C) Response of C BS to decreasing light intensity, under different light 
qualities, estimated by the C 4 model. Error bars represent one SE (n = 6). 



ATP cost of assimilation 

Two empirical ATP costs of (net and gross) assimilation were 
derived. Figure 4 shows the empirical ATP cost of net assimi- 
lation J ATP IA. This quantity expresses the ATP cost involved 
in the assimilation of C0 2 ; that is, how much ATP the plant 
has to produce (= consume, at steady state) to assimilate one 
C0 2 molecule. Figure 4 shows clearly that J ATP IA for HLLL 
plants was very similar to that of LL plants and significantly 
lower than that of HL plants. This means that re-acclima- 
tion was extremely effective in reducing J ATP IA, particularly 
in the vicinity of the LCP Figure 5 shows the ATP cost of 
gross assimilation, J ATP IGA. This quantity is the biochemi- 
cal conversion efficiency of C 4 assimilation, or how much 
ATP is needed to convert bicarbonate into stable assimilates. 
The empirical values for J ATP IGA (Fig. 5, symbols in panels 
A-C) were close to 5.4 and not significantly influenced by 
light intensity or by the growth light regime. This means that, 
in contrast to the common interpretation, the biochemical 
conversion efficiency was not affected by instantaneous light 
intensity. 

To support this result theoretically, we predicted ATP/G^4 
with two different approaches. These methods are compared in 
Fig. 5. A simplified method used 0 as a sole proxy for C 4 oper- 
ating efficiency (cp approach; Fig. 5, solid squares), whereas the 
complete biochemical method (B approach; Fig. 5, solid circles) 
summed the individual ATP demands of processes involved 
in assimilation. Under low light intensities the cp approach 
resulted in an overestimation of J ATP IGA, especially in HL 
plants (Fig. 5A), which display the characteristic hyperbolic cp 
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Fig. 3. Leakiness (<f>) resolved from online isotopic discrimination during 
photosynthesis (A) by means of a full isotopic discrimination model for HL 
plants (squares), LL plants (triangles), and HLLL plants (diamonds). Error 
bars represent one SE (n = 6). 



Table 2. Physiological responses for plants grown under high light (HL), low light (LL) or LL following HL (HLLL) 

The LCP was determined by fitting light curves with dedicated software; Flight was determined by linear regression of A against PAR-Y(ll)/3; BS 
conductance (g BS ) was determined by fitting a modelled J M0D to the measured J ATP (Fig. 3). Different letters identify significant differences across 
rows at P < 0.05 in a Tukey multiple comparison test (Genstat). Mean values ± SE are shown; n = 6 per treatment. 



Unit Mean HL LL HLLL 

LCP |iE-m- 2 -s- 1 15 - 3 24.4 + 1.9 a 10.4 + 0.65 b 11.2 + 1.0 b 

Rught |imol-rrr 2 -s- 1 °- 680 1.05 + 0.14 a 0.510 + 0.057 b 0.477 + 0.053 b 

g ss mol-m- 2 -s- 1 0.000944 0.00136 + 5.2 x 1Q- 4a 0.000647 + 9.2 x 1Q- 5a 0.000822 + 1.9 x1Q- 4a 
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Fig. 4. Measured ATP cost of net assimilation (J AT p/A) for HL plants 
(squares), LL plants (triangles), and HLLL plants (diamonds). Error bars 
represent one SE [n = 6). 
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Fig. 5. ATP cost of gross assimilation, representing C 4 biochemical 
operating efficiency for HL plants (A), LL plants (B), and HLLL plants (C). 
The empirical values for J ATP /GA (empty symbols) were compared to 
predicted values for ATP/G/A (solid symbols) calculated with two different 
approaches. ATP/GA was calculated using <f> as the sole proxy for 
operating efficiency (<j> approach; solid squares) or using a comprehensive 
calculation summing the ATP cost of all processes contributing to 
assimilation (B approach; solid circles). Note that both calculations were 
based on the same dataset, presented in Figs 1-3 and Table 2. Error bars 
represent one SE; n = 6 plants per condition. 



increase in proximity of the compensation point (Fig. 3). Under 
higher irradiances or when (p was lower, the (p approach resulted 
in an accurate estimation of J ATP /GA. The B approach provided 
a better estimate of J ATP /GA, across the range of incident light 
intensities and independent of the values for (p. It is worth stress- 
ing that both these estimates were based on the same dataset 
shown in Fig. 3, but, although the (p approach translated the cp 
pattern shown in Fig. 3 directly into ATP cost, the B approach 
considered the rates of the underpinning biochemical reactions 
and summed the ATP costs involved in each individual process. 

Discussion 

Maize plants were grown under high light then re-acclimated 
to diffuse low irradiance, and compared to plants grown either 
under high light or low light. The particular conditions of re- 
acclimation were intended to represent the transition from full 
sunlight to shaded conditions that maize leaves undergo when 
overgrown by newly emerging leaves at the top of the canopy. 
This is a natural acclimation process for maize leaves, in fact 
for the experimental plants all leaves grown under HL were 
retained throughout the re-acclimation period, and continued 
to photosynthesize under low irradiance, although we cannot 
exclude the possibility that such acclimation to low light could 
also include changes in development or source/sink relations. 
The opposite acclimation is not likely to be a physiologically 
realistic process, and, in fact, when LL plants were moved to 
higher light intensities they promptly shed all leaves grown 
under low light. The natural re-acclimation to low light brought 
about substantial physiological changes, which have implica- 
tions for the energy balance of leaves within a growing canopy. 

Acclimation strategies 

The three types of plants were subject to concurrent gas 
exchange, variable fluorescence, and isotopic discrimination 
measurements. Direct estimates for J ATF were derived from 
a combined low 0 2 -ETR method, and leakiness was derived 
from isotopic discrimination. This comprehensive ecophysio- 
logical characterization highlighted two main re-acclimation 
strategies. 

A first strategy involved reducing R L1GNT . This reduction is 
underpinned by considerable changes at the metabolic level 
that result in reducing the level of basal metabolism. This 
had a direct effect on the LCP and, for this reason, it had a 
direct effect on the ATP cost of net assimilation (see below). 
A second strategy involved the reduction of leakiness (cp). 
HLLL plants showed reduced values for (p as compared to 
HL plants. However the <p hyperbolic increase under low irra- 
diance was similar to that of HL plants. In contrast, and in 
agreement with recent results (Bellasio and Griffiths, 2013; 
Ubierna et al, 2013), LL plants showed a linear trend, with 
(p values that were only marginally affected by irradiance. 
Further model outputs suggested that the general reduction 
in (p observed for HLLL plants (Fig. 3) could only be partially 
accounted by changes in g BS (Table 2), and it may therefore 
be regulated at the level of relative rates of V P and V c , as we 
have discussed previously (Bellasio and Griffiths, 2013). The 



observed plasticity in (p may then involve tuning of biochemi- 
cal reaction rates and, in particular, the ratio between the 
CCM activity and the C 3 activity, or the capacity to accom- 
modate g BS in response to light intensity, as we have recently 
hypothesized (Bellasio and Griffiths, 2013). However, the ulti- 
mate nature of this tuning is still speculative. 

Overall, these strategies were highly effective in reducing the 
ATP cost of net assimilation J ATP IA. In fact, under a PAR of 
50 uE-m~ 2 -s _1 JatpIA for HLLL plants was 35% lower than 
that of HL plants and very similar to that of LL plants (Fig. 4). 
However, this ATP cost reduction was largely associated with 
the reduced R UGH t- In fact, when the effect of Rlight reduc- 
tion was isolated and the biochemical operating efficiency (i.e. 
the ATP cost of gross assimilation J ATP /GA) was considered, 
only minor energetic differences could be observed between 
different light treatments. Re-acclimation significantly influ- 
enced neither the empirical J ATf /GA (mean of 5.47 for HLLL 
and 5.45 for HL, as compared to 5.40 for LL) nor the predicted 
ATP cost of GA (ATP/GA, B approach; mean 5.56 for HLLL 
and 5.42 for HL, as compared to 5.61 for LL). This shows that 
if there were any effect of varied <p on the overall biochemi- 
cal conversion efficiency the effect was undetectable using the 
methods described. On one hand this confirms the difficulties in 
estimating leakiness from leaf-level energetics (Furbank et al. , 
1990; Kromdijk, 2010), on the other it highlights the complex- 
ity of the leakiness phenomenon, which depends at the same 
time on anatomical and biochemical traits. In this study we 
have specifically addressed the ATP demand, but other aspects 
are intertwined and may all contribute to (p dynamics. These 
could include (von Caemmerer and Furbank, 2003; Furbank, 
201 1; Bellasio and Griffiths, 2014) (i) regulating the ratio of C 4 
dicarboxylic acid to amino acid export to BS, (ii) regulating 
reducing power export from mesophyll to BS cells in response 
to demand, (iii) partitioning metabolic work between contrast- 
ing cells types (e.g. PGA reduction, starch synthesis, glycolate 
recycling, RuBP + PEP regeneration), (iv) optimizing energy 
availability in BS and mesophyll cells while at the same time 
(v) maintaining the equilibrium between the CCM and the C 3 
activity, and, finally, (vi) trading-off at the level of BS conduct- 
ance, between the capacity to support very high diffusion (and 
assimilation) rates and the necessity to limit leakage of C0 2 
out of the BS (Sowinski et al , 2008). 

Predicting C 4 operating efficiency 

The 'conventional' approach to predicting C 4 biochemical 
operating efficiency (i.e. the ATP cost of gross assimilation, 
ATP/GA) uses leakiness, cp as the sole proxy (eqn 2, referred 
as the 4> approach). With the (p approach the C 3 activity is 
considered to have an invariable cost of 3 ATP/GA (pho- 
torespiration is neglected), whereas the CCM is assumed to 
be supplied solely by PEPC activity, which is assumed to refix 
entirely the retrodiffused (leaked) C0 2 . Our empirical evi- 
dence largely confirms the validity of the cp approach, which 
closely predicted the trend and the magnitude of J ATP /GA 
under PAR > 125 uE-m~ 2 -s _1 for LL + HLLL plants, and 
under PAR >500 uE-m~ 2 -s _1 for HL plants. However, under 
low irradiances, the cp approach overestimated the trend of 
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J ATP IGA, especially for HL plants. This overestimation is 
dependent on the assumptions of the cp approach (photores- 
piration is neglected and the CCM is driven solely by PEPC), 
which hold only under high irradiances, while under low irra- 
diance they are no longer valid. In fact, PEPC and Rubisco 
activities proportionally decrease under decreasing irradi- 
ance, limited by the decreasing ATP availability. As opposed 
to that, BS respiration is largely unaffected by light intensity, 
and, under decreasing light intensities, the BS-respired C0 2 
progressively outweighs PEP carboxylation rate (V P ). Hence, 
cp — that is, the ratio of retrodiffusing C0 2 over PEP carboxy- 
lation rate — becomes progressively higher as light intensity 
approaches the compensation point. This gives rise to the 
extensively documented (empirically and theoretically) hyper- 
bolic (p increase (for review see Ubierna et al, 2011), which 
can be largely supplied by respiration without an additional 
engagement of PEPC. A constant degree of engagement of 
PEPC, even under the hyperbolic cp increase, is consistent 
with the observation that both the ratio of PEPC/Rubisco 
carboxylation rate (VplV c ) and the optimal partitioning 
factor between the CCM activity and the C 3 activity (x) are 
largely independent of light intensity (von Caemmerer, 2000; 
Kromdijk et al., 2010). 

These considerations can be better appreciated if the CCM 
is viewed as complex machinery. The activity of PEPC is 
only one of the systems which contribute to loading C0 2 
into the BS. Recently it has become increasingly clear how 
photorespiration may contribute to the CCM, and is the pre- 
dominant driving force in evolutionally early types of CCM 
(Sage et al, 2012; Schulze et al, 2013). We have shown in this 
paper and in previous work (Bellasio and Griffiths, 2013) that 
the CCM can be increasingly supplied by respiration under 
limiting light conditions, bringing about increased leakiness 
even without a predicted increase in the activity of PEPC. 
It is worth remembering that the compartmentalization of 
photochemical water oxidation to mesophyll cells, whose 
degree may vary considerably between subtypes and along 
the evolutionary line (Meierhoff and Westhoff, 1993; Sage, 
2004; Furbank, 2011; Sage et al, 2011), also contributes to 
increasing the ratio of 0 2 /C0 2 at the active site of Rubisco 
and should also be considered as a component of the com- 
plex machinery of the CCM. 

In view of this complexity, leakiness, which reflects inher- 
ently complex biochemical and anatomical traits, should only 
be used to predict the magnitude (and ATP cost) of the CCM 
under high light regimes. However, we showed that J ATP IGA 
could be closely predicted using a complete biochemical 
approach (B approach, eqn 2), whereby the ATP cost of 
all processes contributing to assimilation are summed. We 
used the equations that we have recently derived (Bellasio 
and Griffiths, 2014), which are based on the comprehensive 
description of the C 4 metabolism outlined by Furbank (201 1) 
and on the validated C 4 model (von Caemmerer, 2000, 2013). 
Within this approach, the ATP-consuming processes consid- 
ered are PGA reduction, PEP regeneration (through PPDK 
and PEPCK), and starch synthesis. Furthermore the B 
approach subtracts the PGA used by respiration, which does 
not need to be reduced (PGA reduction to DHAP consumes 
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1 ATP and 1 NADPH). With this comprehensive calculation, 
the empirical data could be closely predicted in the vicinity of 
the LCP. Notably, under decreasing light intensities, the bio- 
chemical conversion efficiency did not decrease, regardless of 
the hyperbolic (p increase observed in HL and HLLL plants. 



Conclusion 

In this study we set out to investigate the strategies deployed by 
maize plants grown under high light intensities when re-accli- 
mated to low light. We showed that the main re-acclimation 
drivers were the reduction of respiration and the reduction of 
leakiness, and these were likely to be accompanied by complex 
metabolic reorganization. Overall, these strategies were very 
effective in reducing the ATP cost of net assimilation under 
low light intensities, which, for HLLL plants, decreased by 
35% as compared to HL plants under PAR = 50 u.E-m~ 2 -s _1 . 
This shows clearly that the net energy conversion efficiency 
under limiting light is to a considerable extent ameliorated by 
the acclimation of mature leaves to low light. 

By calculating ATP cost of gross assimilation we could iso- 
late the contribution of day respiration from the other bio- 
chemical effects (which include the reduction of leakiness). The 
ATP cost of gross assimilation was not significantly different 
for HLLL plants as compared to HL plants. This showed that 
re-acclimation did not change the efficiency of C 4 metabolism, 
even if it considerably reduced leakiness, implying that the 
effect of reduced leakiness on C 4 energetics was not detectable. 
Leakiness dynamics may then be associated to other processes 
occurring at biochemical level such as the regulation between 
BS versus mesophyll metabolic engagement and CCM versus 
C 3 activity. In addition, we provided compelling theoretical and 
empirical evidence showing that the increase in hyperbolic leak- 
iness, observed under low light intensities (Fig. 3), is not associ- 
ated with a loss of energetic efficiency. The well-consolidated 
idea of C 4 efficiency loss under low light conditions (e.g. von 
Caemmerer, 2000; Tazoe et al, 2008) relies on assumptions that 
should be reconsidered in view of recent discoveries: the CCM 
is not uniquely supplied by the ATP-costly PEPC activity, but, 
under certain conditions, the contribution through respiration 
and photorespiration may be significant (Kromdijk et al. , 2010; 
Sage et al, 2012; Bellasio and Griffiths, 2013). We proposed 
a comprehensive biochemical method (Bellasio and Griffiths, 
2014) based on the validated C 4 model (von Caemmerer, 2000). 
The biochemical method predicts the C 4 conversion efficiency 
(as ATP cost of gross assimilation), taking into account the 
active and passive contributions to the CCM. 

The implications for loss of productivity at the field 
scale being specifically associated with increased leakiness 
(Kromdijk et al, 2008) may be less severe than previously 
thought. However, here we have shown the potential for accli- 
mation with a somewhat extreme acclimation pattern whereby 
mature leaves were switched from the full light to deep shade. 
Realistically, mature leaves will undergo a more gradual tran- 
sition from full sunlight through a condition characterized 
by rapid changes in irradiance (daily shading, sunflecks), to 
complete shade. The actual extent to which leaves optimize 



energy efficiency, when exposed to such a complex pattern of 
illumination under field conditions, remains to be addressed. 

Supplementary material 

Supplementary material is available at JXB online. 

Supplementary Fig. SI. | values for the calculation of A 
(equations are reported in Table 1). 
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